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Abstract 

Azimuthal asymmetries in exclusive electroproduction of a real photon from a longitudinally 
polarized deuterium target are measured with respect to target polarization alone and with 
respect to target polarization combined with beam helicity and/or beam charge. The asym- 
metries appear in the distribution of the real photons in the azimuthal angle <j> around the 
virtual photon direction, relative to the lepton scattering plane. The asymmetries arise from 
the deeply virtual Compton scattering process and its interference with the Bethe-Heitler pro- 
cess. The results for the beam-charge and beam-helicity asymmetries from a tensor polarized 
deuterium target with vanishing vector polarization are shown to be compatible with those 
from an unpolarized deuterium target, which is expected for incoherent scattering dominant at 
larger momentum transfer. Furthermore, the results for the single target-spin asymmetry and 
for the double-spin asymmetry are found to be compatible with the corresponding asymmetries 
previously measured on a hydrogen target. For coherent scattering on the deuteron at small 
momentum transfer to the target, these findings imply that the tensor contribution to the cross 
section is small. Furthermore, the tensor asymmetry is found to be compatible with zero. 
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1. Introduction 

Generalized Parton Distributions (GPDs) provide a framework for describing the mul- 
tidimensional structure of the nucleon [1-3]. GPDs encompass parton distribution func- 
tions and elastic nucleon form factors as limiting cases and moments, respectively. Parton 
distribution functions are distributions in longitudinal momentum fraction of partons 
in the nucleon, and are extracted from measurements of inclusive and semi-inclusive 
deep-inelastic scattering. Form factors are related to the transverse spatial distribu- 
tion of charge and magnetization in the nucleon. Both form factors and (transverse- 
momentum-integrated) parton distribution functions represent one-dimensional distri- 
butions, whereas GPDs provide correlated information on transverse spatial and longitu- 
dinal momentum distributions of partons [4-9]. Furthermore, access to the total parton 
angular momentum contribution to the nucleon spin may be provided by GPDs through 
the Ji relation [3]. 

Hard exclusive leptoproduction of a meson or photon, with only an intact nucleon or 
nucleus remaining in the final state, can be described in terms of GPDs. GPDs depend on 
four kinematic variables: t, x, £, and Q 2 . In this case, t is the Mandclstam variable, or the 
squared four-momentum transfer to the target, given by t = (p — p') 2 , where p (p 1 ) is the 
initial (final) four-momentum of the target. In the 'infinite' target-momentum frame, x 
and £ are related to the longitudinal momentum of the parton involved in the interaction 
as a fraction of the target momentum. The variable x is the average momentum fraction 
and the variable £, known as the skewness, is half the difference between the initial and 
final momentum fractions carried by the parton. The evolution of GPDs with Q 2 = — q 2 , 
with q = k — k! the difference between the four-momenta of the incident and scattered 
leptons, can be calculated in the context of perturbative quantum chromodynamics as in 
the case of parton distribution functions. This evolution has been evaluated to leading 
order [1-3,10] and next-to-leading order [11-13] in the strong coupling constant a s . The 
skewness £ can be related to the Bjorken scaling variable xb = Q 2 /(2p ■ q) through 
£ ~ xb/(2 — xb) in the generalized Bjorken limit of large Q 2 , and fixed ig and t. There 
is currently no consensus as to how to define £ in terms of experimental observablcs; 
hence the experimental results arc typically reported as projections in xb- The entire 
x dependences of GPDs are generally not experimentally accessible, an exception being 
the trajectory x = £ [14,15]. 

GPDs can be constrained through measurements of cross sections and asymmetries 
in exclusive processes such as exclusive photon or meson production. In this paper, the 
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Deeply Virtual Compton Scattering (DVCS) process, i.e., the hard exclusive production 
of a real photon, is investigated using a longitudinally polarized deuterium target. 

The spin-1/2 nucleon is described by four leading-twist quark-chirality-conserving 
GPDs H, E, H and E [1-3,16]. In contrast, DVCS leaving the spin-1 deuteron intact re- 
quires nine GPDs: Hi, H 2 , H 3 , H4, H 5 , Hi, H 2 , H 3 and H4 [17-19]. In the forward limit 
of vanishing four- momentum transfer to the target nucleon (t — > and £ — > 0), the pairs 
of GPDs (H, Hi) and (H, Hi) reduce respectively to quark number density and helicity 
distributions. In this limit the GPD H$, sensitive to tensor effects in the deuteron, re- 
duces to the tensor structure function bi , which was measured in inclusive deep-inelastic 
scattering on a tensor polarized deuterium target [20]. Both H3 and H$ are associated 
with the 5% D-wave component of the deuteron wave function in terms of nuclcons [22] . 
In addition to GPD Hi, they both contribute to the beam-helicity and beam-charge 
asymmetries. The term with GPD H§ dominates in the beam-hclicity<53tensor asymme- 
try in DVCS from a longitudinally polarized deuterium target at very small values of 
t [18]. At this kinematic condition, the asymmetry with respect to target polarization is 
dominated by the term with GPD Hi. Thus, the measurement of certain asymmetries in 
DVCS on a polarized deuterium target may provide new constraints for these GPDs. 

This paper reports the first observation of azimuthal asymmetries with respect to 
target polarization alone and with respect to target polarization combined with beam 
helicity and/or beam charge, for exclusive electroproduction of real photons from a lon- 
gitudinally polarized deuterium target. The asymmetries arise from the DVCS process 
where the photon is radiated by the struck quark, and its interference with the Bethe- 
Heitler (BH) process where the photon is radiated by the initial or final state lepton. 
The resulting asymmetries combine contributions from the coherent process ed — > ed^y, 
and the incoherent process e d — > epnj where in addition a nucleon may be excited to a 
resonance. The coherent reaction contributes mainly at very small values of t, while the 
incoherent process dominates elsewhere. It is natural to model the incoherent process as 
scattering on only one nucleon in the deuteron, while the other nucleon acts as a spec- 
tator. Monte Carlo simulations in HERMES kinematic conditions [23] suggest that the 
proton contributes about 75% of the incoherent yield and the neutron about 25%, and 
included in these, nucleon resonance production contributes about 22% of the incoher- 
ent yield. The incoherent reaction on a proton dominates that on a neutron because of 
the suppression of the BH amplitude on the neutron by the small elastic electric form 
factor at low and moderate values of the momentum transfer to the target. The depen- 
dence of the measured asymmetries on the kinematic conditions of the reaction is also 
presented and these results on the deuteron are compared where appropriate with the 
corresponding results obtained on a longitudinally polarized hydrogen target [24] . 

2. Deeply virtual Compton scattering 

2.1. Scattering amplitudes 

The DVCS process is currently the simplest experimentally accessible process that can 
be used to constrain GPDs. The initial and final states of DVCS are indistinguishable 
from those of the competing BH process. For a target of atomic mass number A and 
no target polarization component transverse to the direction of the virtual photon, the 
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general expression for the cross section of the coherent reaction ei-> e A 7 or incoherent 
reaction ei->e(4 - l)N~f reads [18,25] 

da x A e 6 |T| 2 



cLm dQ 2 d|*| d</> 32 (2tt) 4 Q 4 ^ ' ^ 

Here, xa = Q 2 / '(2M a v ) is the nuclear Bjorken variable, where Ma is the mass of the 
nucleus and v = p ■ q/MA, e is the elementary charge, e = 2xaMa/ -y/Q 2 " and T is the 
reaction amplitude. The azimuthal angle of the real photon around the virtual-photon 
direction, relative to the lepton scattering plane, is denoted by 4>. The cross section 
contains the coherent superposition of BH and DVCS amplitudes: 

\T\ 2 = |7bh + 7dvcs| 2 = |7bh| 2 + |7dvcs| 2 + Tdvcs 7bh + Tdvcs ^bh , (2) 

" v ' 

I 

where T denotes the BH-DVCS interference term. The BH amplitude is calculable to 
leading order in QED using the form factors measured in elastic scattering. 

The interference term I and the squared DVCS amplitude |7dvcs| 2 in Eq. 2 provide 
experimental access to the (complex) DVCS amplitude through measurements of various 
cross section asymmetries as functions of <j> [18]. Each of the three terms of Eq. 2 can 
be written as a Fourier series in </>. In the case that the beam and the target may be 
longitudinally polarized, these terms read 

|T ™ l ^^vltw !< S c " BHc " ( "* , ■ (3) 



n=0 



IWsl 2 = AWs x i ^ VCS cos(n0) + J2 *n VCS M™t>)\ 

lri=0 n=l J 



(4) 



The symbols A'bh = ■> ttt 2^ 1 A'dvcs = tti and K\ = — denote kinematic factors, 
where y = p ■ q/(p ■ k), and eg stands for the (signed) lepton charge in units of the ele- 
mentary charge. In the case of unpolarized beam and target, certain coefficients vanish. 
All Fourier coefficients c n and s n in Eqs. 3-5 depend on the longitudinal target polar- 
ization, with some also having a dependence on the beam helicity. The coefficients c^ H 
in Eq. 3 depend on electromagnetic form factors of the target, while the DVCS (interfer- 
ence) coefficients c^ vcs (cjj and s^ vcs (sjj involve various GPDs. The squared BH and 
interference terms in Eqs. 3 and 5 have an additional <j> dependence in the denominator 
due to the lepton propagators Pi{4>) and Vi(<$) [25,16]. The Fourier coefficients cj> and 
in Eq. 5 can be expressed as linear combinations of Compton Form Factors (CFFs) [18], 
while the coefficients cj^ vcs and s^ vcs arc bilinear in the CFFs. Such CFFs are convo- 
lutions of the corresponding GPDs with the hard scattering coefficient functions. 

For a longitudinally (L) polarized lepton beam scattered from an unpolarized tar- 
get, the beam-charge asymmetry Ac and the charge- difference beam- helicity asymmetry 
A\aj (sensitive to the interference term) and charge-average beam-helicity asymmetry 
-^lu CS (sensitive to the squared DVCS term) can be measured if all four combinations 
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of beam charge and helicity are available [26] . Results of their Fourier amplitudes for an 
unpolarized deuterium target were recently published by HERMES [27]. 

Unfortunately, the present data set for a longitudinally polarized target does not in- 
clude all four combinations of beam charge and sign of beam polarization. Therefore, the 
beam-helicity asymmetries presented in this paper are single- charge observables, which 
entangle the interference and squared DVCS term. Fortunately, measurements of charge- 
averaged beam-helicity asymmetries on hydrogen [26] and deuterium [27] targets showed 
that the contribution by the squared DVCS term is negligible in HERMES kinematic 
conditions, at the precision of these measurements. 



2.2. DVCS on the deuteron 

For coherent scattering on a spin-1 target nucleus polarized longitudinally with re- 
spect to the virtual photon direction, and with spin projection A = ±1,0, the following 
decomposition of the Fourier coefficients appearing in Eqs. 3-5 is introduced [18]: 



c*(A) = ^A 2 c «; unp + Ac£ LP + (1 - ^A 2 )c«; LLP 



with R e {BH,DVCS,I}, and similarly for the s„ coefficients with R e {DVCS, I}. The 
subscript 'unp' denotes unpolarized and 'LP' and 'LLP' denote respectively vector and 
tensor terms for parts of the cross section related to longitudinal polarization. For an 
unpolarized target nucleus, one recovers the value [c^ (A = —1) + c^(A = 0) + cJJ(A = 
+l)]/3 = c^ unp . A purely tensor-polarized target nucleus with A = results in = 
c nLLP' wm le for A ^ all coefficients contribute. 

Equation 6 is applicable only for purely polarized states with A = ±1,0. In a real 
experiment, the longitudinally polarized deuterium target contains a mixture of these 
pure polarized states, characterized by vector and tensor polarizations P z and P zz defined 



n + — n n + + n — 2n° 

2 = n+ + n-+n° ' zz = n+ + n~ + n° ' ™ 



where n + , n~ and n° arc the populations of the state with A = +1, —1 and 0, respectively. 

For a lepton beam with given longitudinal beam polarization Pi scattering coherently 
on a deuterium target with given vector and tensor polarizations P z and P zz , the Fourier 
series of the squared reaction amplitude reads, using the spin decompositions of Eq. 6, 
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|7bh| 2 



K 



BH 



i 



BH 



,LLP 



) cos(n^) 



2 zz 



2 



E C «,unp C Os(n0) 
n=0 

1 2 

^ E C «-LP C0S M) + 2 P " E( C n"np - 

n=0 n=0 

Tdvcs I 2 = #dvcs I E Cp cos (^) + p z s ?,u^p sin 

I n=0 

^E c «X? S cos(^) + ^^vcs sin(r 

1=0 n=l 

E(C? P S - cos(n</>) + ft ( S ?vc p s - s ?vcs } sin 

( 3 2 

E c Un P cos (™0) + p! E s «m P sin(n0) 

n=0 n=l 
3 

p t E C «>LP cos(n0) + Y s -i,lp sin(n0) 

i=0 n=l 
3 2 

E( C «^"P _ C »,LLP) COS (^) + P i ^i^unp ~ 4,LLp) sin(n</>) 



P* 



n=0 



(8) 



(9) 



(10) 



Note that the beam polarization Pj? and the target vector and tensor polarizations P z 
and P zz are here factored out of the corresponding Fourier coefficients in Eqs. 3-5, thus 
leaving only the dynamical kinematic dependences encoded in the Fourier coefficients in 
Eqs. 8-10. 



2.3. Asymmetries on the deuteron 

For data with longitudinal polarization of both beam and target, the following nota- 
tion is introduced: — > (■<— ) to denote positive (negative) beam helicity, and =>■ and <= 
to denote the deuteron target vector-polarization direction anti-parallel and parallel to 
the beam momentum direction in the target rest frame. In contrast to lepton scattering 
off longitudinally polarized hydrogen [24], there are many more observables (asymme- 
tries) in the case of deuterium. They may be classified according to whether the cross 
section for A = explicitly appears in the definition of this asymmetry. An example of 
the 'incomplete' asymmetries where it does not appear is the beam-helicity asymmetry 
A^(e£, P zz , (j)): defined for beam charge eg and tensor polarization P zz as 



A(j^{ei,P zz ,4>) 


+ da^{ee.,P zz ,4>) 




da^{e e ,P zz ,0) 


+ da^{e e ,P ZZl (j>) 


da^(e e ,P zz ,(/}) 


+ d<r^(e e ,P zz ,(j)) 


+ 


da^(e,,P Z2 ,0) 


+ da^(e e ,P zz ,(j)) 



(11) 



Here, the symbol 'do - ' denotes a generic differential cross section. 
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For coherent scattering, and to leading order in a s and in leading twist, the expansion 
in powers of the Bjorken variable xu for the deuteron target, and r = i/(4M|,), where 
Me is the deuteron mass, yields [18] 



-i. 



x D {2-v)^{l-y) 



2-2y + y 2 
G\Hi — ^GiH§ — t G\H 3 + G 3 (T-Li — 5^5) 



sun 



2t 2 G 3 H 3 



G 2 -2tG 1 G 3 + 2t 2 G 2 



-<-'( 



iy - 



sin < 



(12) 
(13) 



Here, G\ and G3 are deuteron elastic form factors [28]. (For comparison with experimental 
data, the actual value of P zz 7^ 1 must be taken into account in, e.g., Eqs. 12, 13, 
and those that follow.) Equation 13 is obtained neglecting the contributions of non- 
leading terms in r in Eq. 12, which are less than 10% at —t < 0.03 GcV 2 (see Fig. 3 in 
Ref. [27]). As can be seen from Eqs. 12 and 13, this asymmetry involves a different linear 
combination of the imaginary parts of the deuteron CFFs Hi, H 3 and H5 compared to 
the asymmetry ^1^(0) (see Eqs. 25-27 in Ref. [27]). More specifically, any difference 
between these two asymmetries at small values of —t may be ascribed to the CFF T-L^. 
Detailed information about the relations between these CFFs and corresponding GPDs 
can be found in Ref. [18]. 

Similarly, the beam-charge asymmetry for tensor polarization P zz is defined as 





+ &(T t {P zz ,4>) 




da^(P zz ,<t>) 


+ da < 


-{Pzz A) 


da^(P zz ,4>) 


+ da^{P zz ,4>) 


+ 


da^{P zz ,4») 


+ der< 


-{Pzz A) 



, (14) 



where the symbols + (— ) denote positive (negative) beam charge. For coherent scattering, 
the cos 4> component in the kinematic expansion of Eq. 14 is sensitive to the real part of 
the same linear combination of CFFs as that appearing in Eq. 13: 



A c J=(Pzz = +!•< 



XD 



y 



■ cos ( 



(15) 



The different sign of the asymmetry A C ^(P ZZ = +1,0) compared to Ref. [18] is due to 
the use of the Trcnto convention [29] in this work, i.e., <fi = n — 0hgi- 

Another single-charge beam-helicity asymmetry, which differs from A^s. (eg, P zz , cj>) 
and A c <^(eg, P zz , (f>), involves polarized beam and (longitudinal) tensor polarization of 
the deuteron: 

A Lzz (eg, 4>) = - - - , , - , , , , (16) 



3der~ 
2d(7° and da, 



i P (e£» + 3dcr^ p (e f » ' 

■|(dcr^ +d(T <;= + der°), where da represents 



with dfi zz = der^ +da *, uu ^ unp — 3 

the cross section for deuterons in the A = state. For coherent scattering, the asymmetry 
Ahzz(ee, 4 1 ) involves a different linear combination of the imaginary parts of the deuteron 
CFFs Hi, H 3 and H 5 compared to A^(eg, P zz = +1, 4>) and A l w (<j)): 
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v4-Lzz(e 



2x D (2-y)yJ$(l-y) 



sin ( 



2-2y + y 2 

dHs + t (Gi^ 3 + - § G 3 H 5 ) - 2t 2 G 3 U 3 
1 3G 2 -4rGiG 3 + 4t 2 G 2 



'■ ('1 



2x D (2~y)^(l~y)c s . 



3Gi 



■ sine 



(17) 
(18) 



2 - 2y + y 2 

Finally, the single-charge asymmetry with respect to longitudinal vector polarization 
of the target is defined as 



A\jh(e£, P Z z,4>) = 



da^(e e ,P zz ,0) 


+ d<j^(e e ,P zz ,(j)) 




da^(e£,P zz ,4 


)+da^(e t ,P zz ,(j)) 


da^(e e ,P zz ,^) 


+ da^{e e ,P zz ,</)) 


+ 


da^(e e ,P zz ,4 


) + d<T^(ei,P zz ,(t>) 



(19) 



For coherent scattering, in analogy to the previously elaborated asymmetries, it reduces 
to 



A\jh(ee, Pzz = +1, 4>) — ~ e e- 



x D 



y 



x inn 



G\H.\ + G2 {Hi — \Hz 



t I G 3 T~L\ + ^p- G-{H.^ 



G'f - 2rGiG 3 + 2r 2 G 2 



xd Vg(l - y) 9m G{H 1 + ^G 2 (Hi - 



(20) 
(21) 



y Gj 
Thus, this asymmetry is sensitive to the imaginary part of the Compton form factor Hi 



3. The HERMES experiment 



A detailed description of the HERMES spectrometer can be found in Ref. [30] . A lon- 
gitudinally polarized positron or electron beam of energy 27.6 GeV was scattered off a 
longitudinally polarized deuterium gas target internal to the HERA lepton storage ring 
at DESY. The lepton beam was transversely self-polarized by the emission of synchrotron 
radiation [31]. Longitudinal polarization of the beam at the target was achieved by a pair 
of spin rotators in front of and behind the experiment [32]. The sign of the beam po- 
larization was reversed approximately every two months. Two Compton backscattering 
polarimeters [33,34] measured independently the longitudinal and transverse beam po- 
larizations. The average values of the beam polarization for the various running periods 
are given in Table 1; their average fractional systematic uncertainty is 2.2%. 

The target cell was filled with nuclear-polarized atoms from an atomic beam source 
based on Stern Gcrlach separation with radio- frequency hypcrfinc transitions [35]. The 
polarization and atomic fraction of the target gas were continuously monitored [36,37]. 
Most of the longitudinally polarized deuterium data were recorded with average vector 
polarizations 0.851 ± 0.031 and —0.840 ± 0.028, and with an average tensor polarization 



9 



Table 1 

The sign of the beam charge, the luminosity-averaged beam polarization and target vector and tensor 
polarization values, for the years 1998-2000 and the integrated luminosity of the data sets used for the 
extraction of the various asymmetry amplitudes (see Table 3) on a longitudinally polarized deuterium 
target. The uncertainties for the polarizations are given in the text. 



Lcpton 


Beam 


Target 




Luminosity [pb 1 


] 




Year Charge 


Polarization 


Polarization 


















P z 


Pzz 


(Ail, ^ll) 


{A ,A C 


) 

L 




1998 e- 


- 0.509 




±0.856 


f 0.827 




26.2 






1999 e+ 


- 0.547, ± 


0.518 


±0.832 


f 0.827 


29.7 


14.2 




29.7 


2000 e+ 


-0.537, ± 


0.524 


-0.840, ±0.851 - 


f 0.827 


125.8 


43.5 




125.8 


2000 e+ 


-0.542, + 


0.525 


-0.010 


- 1.656 








22.7 


Sum 










155.5 


83.9 




178.2 



Table 2 

The integrated luminosity of the data used for the extraction of various asymmetry amplitudes (see 
Table 3) on a longitudinally polarized deuterium target for each lepton beam charge and sign of the 
polarization. 



Lcpton Charge 


Sign of the Beam Polarization 


Luminosity [pb x ] 


e~ 


negative 


26.2 


e~ 


positive 




e+ 


negative 


75.4 


e+ 


positive 


102.8 



of 0.827 ± 0.027 [38] (corresponding to a small population of the A = state). The 
extraction of Ahzz{4>) employed the fraction of the data taken in the year 2000 recorded 
with a tensor-polarized deuterium target where deuterons in the A = state were injected 
into the target cell, resulting in an average tensor polarization of —1.656 ± 0.049 with 
negligible vector polarization (— 0.010±0.026). The amount of data accumulated for each 
lcpton beam charge and sign of the polarization are summarized in Table 2. 

The scattered leptons and produced particles were detected by the HERMES spectrom- 
eter in the polar angle range 0.04 rad < 9 < 0.22 rad. The average lepton identification 
efficiency was at least 98% with hadron contamination of less than 1%. 

4. Event selection 

The data sets used in the extraction of the various asymmetries reported here are given 
in Table 1. In this analysis, it was required that events contain exactly one charged- 
particle track identified as a lepton with the same charge as the beam lepton, and 
one photon producing an energy deposition E~ > 5 GeV (> 1 MeV) in the calorimeter 
(preshower detector). The following kinematic requirements were imposed on the events, 
as calculated from the four-momenta of the incoming and outgoing lepton: 1 GeV 2 < 
Q 2 < 10 GeV 2 , Wl > 9 GeV 2 , v < 22 GeV and 0.03 < x N < 0.35, where = 
M$f + 2Mnv — Q 2 and xn = Q 2 /(2Mjvf)- For the nucleonic mass Mn, the proton mass 
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was used in all kinematic constraints on event selection, even at small values of —t where 
coherent reactions on the deuteron are dominant, because the experiment does not dis- 
tinguish between coherent and incoherent scattering and the latter dominates over most 
of the kinematic range. Monte Carlo studies have shown that this choice has little effect 
on the result [23]. In order to reduce background from the decay of neutral mesons, the 
angle between the laboratory three-momenta of the real and virtual photons was limited 
to 9-y*-y < 45mrad. The minimum angle requirement # 7 » 7 > 5mrad was chosen accord- 
ing to Monte Carlo simulations in order to ensure that the azimuthal angle <p remains 
well-defined while accounting for the finite angular resolution of the spectrometer. 

An 'exclusive' event sample was selected by requiring the squared missing mass M\ to 
be close to the squared nucleon mass Mjf, where M\ is defined as M\ = (q + P N - q') 2 
with Pjv = (Mjv, 0, 0, 0) and q' the four-momentum of the real photon. The exclusive re- 
gion is defined as -(1.5) 2 GcV 2 < M% < (1.7) 2 GcV 2 to minimize background from deep- 
inelastic scattering fragmentation processes, while maintaining reasonable efficiency [39]. 

As the recoiling target nucleon or nucleus was undetected, the Mandclstam variable 
t was reconstructed from the measured four-momenta of the scattered lepton and the 
detected photon. The resolution in the photon energy from the calorimeter is inadequate 
for a precise determination of t. Hence for events selected in the exclusive region in M\ , 
the reaction is assumed to take place on a nucleon and the final state is assumed to 
contain only the scattered lepton, the real photon and the nucleon that was left intact 
(e N — > eNj). This allows t to be calculated with improved resolution using only the 
photon direction and the lepton four- momentum [40] : 

t = -Q 2 -2v{v- ^ 2 + Q 2 cosfl r7 ) 

1 + ife (" - v 7 ^ 2 + Q 2 cos<v 7 ) 

The error caused by applying this expression to incoherent events with a nucleon excited 
to a resonance in the final state is accounted for in the Monte Carlo simulation that is 
used to calculate the fractional contribution of background processes per kinematic bin. 
This simulation also demonstrated that this method is applicable also to coherent events. 
A further restriction, —t < 0.7 GeV 2 , is used in the selection of exclusive events in order 
to reduce background. 

The exclusive sample comprises coherent and incoherent scattering, including reso- 
nance excitation. Over most of the kinematic range incoherent scattering dominates. 
The events from coherent scattering off the deuteron are concentrated at small values 
of —t. The Monte Carlo simulation showed that requiring —t < 0.06 GeV 2 enhances the 
relative contribution of the coherent process from 20% to 40% in the data sample. Re- 
quiring — t < 0.01 GeV 2 can further enhance the coherent contribution to 66%, but only 
at the cost of a rapidly decreasing yield. The first bin defined in Section 7 covering the 
range —t < 0.06 GeV 2 is sensitive to coherent effects. 

5. Extraction formalism 

The simultaneous extraction of Fourier amplitudes of beam-charge and beam-helicity 
asymmetries combining data collected during various running periods at HERMES for 
both beam charges and hclicities on unpolarizcd hydrogen or deuterium targets is de- 
scribed in Rcfs. [27,26]. It is based on the maximum likelihood technique [41], which 
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provides a bin- free fit in the azimuthal angle (f> (sec Rcf. [42] for details). In this paper, 
data taken with a longitudinally polarized deuterium target were analyzed with a similar 
technique. In the fit, event weights were introduced to account for luminosity imbalances 
with respect to beam charge and polarization. 

Because the target polarization was longitudinal with respect to the direction of the 
incoming beam, the data also contain contributions arising from the small transverse 
polarization with respect to the direction of the virtual photon. This 6% — 12% trans- 
verse component of the target polarization, depending on the kinematic conditions of 
each bin, was neglected in the formalism presented. Hence, the extracted Fourier compo- 
nents contain contributions from this transverse component. However, mainly non-leading 
(higher-twist) amplitudes are affected by this choice. These effects are estimated from 
the measurement of the transverse-target-spin asymmetries at HERMES [42] to be less 
than 0.008 on a proton target, and hence are expected to be negligible compared with 
the uncertainties here. 



5.1. Single- charge formalism 

Data collected with an e~ beam and a polarized deuterium target were not used for the 
extraction of harmonics of A L ^ , A\jl and All because only negative beam polarization 
is available for this charge. Hence, Fourier amplitudes of the three single-charge asym- 



metries A^(ee 



), A VL (e e = +1,P Z2 » and All (fit = + 1,P Z , 



, defined 

respectively in Eqs. 11, 19, and 30, are simultaneously extracted using data from scatter- 
ing of a longitudinally polarized positron beam off a longitudinally polarized deuterium 
target. 

The distribution in the expectation value of the yield can be written as 



d(A0(e £ 



-l,P h P z ,P Zi 

X d(T uS (£ 



C(e e = +l,Pi,P„P„)ri(4>) 



P,A uh (e l = +l,P„,<l>) + P l P z A L L(et = +l,P„,<i>) , (23) 



where C denotes the integrated luminosity and 77 the detection efficiency. The cross 
section for the production of real photons by unpolarized positrons on a tensor-polarized 
deuterium target with vanishing vector polarization is given by 



da v ^(e e = +l,P zz ,<p) = - 



d<7^ + (P zz ,0)+ckr^+(.P zz 



d^ + (P zz ,, 



do* 



-(Pzz,< 



K 



K 



BH 



[E c ™npcos(n0) + -P zz ^(, 



BH 



c n.LLp) COs(n</>) 



n=0 



1 

+ ^DVCS [ £ C^CS cos(n0) + _p zz J2(c 
n=0 

3 



DVCS 
n.unp 



„DVCS 
-n,LLP 



) cos(n</>) 



K 



COs(n(/>) + -P Z z y^(Cn. 



imp 



-n.LLP 



) cos(n^) 



(24) 



(25) 
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where A 



is a common kinematic factor. 



32 (27r) 4 Q 4 v / T+S 7 

The single-charge asymmetries appearing in Eq. 23 are expanded in terms of the same 
Fourier harmonics used in the expansion of the cross section in Eqs. 8-10 and in the 
numerators appearing in Eqs. 26, 28, and 31: 



K 



da uS (e £ = +1,P ZZ ,, 



x |A' D vcs 
K 



,DVCS 
5 l,unp 



1 



cin^+lp f«,DVCS _ DVCS\ • 

sm q> + -f zz (s ljUnp s 1LLP j sin 



2 



imp s n,LLP 



I 



) sm(n(f>) 



s'm(n4 



'{e e = +1, P zz ) sm(nq 



(26) 
(27) 



Ajjh(ei = +1,P ZZ ,4>) = - — — — 

da v ^(e e = +l,P zz ,(p) 

x | A'dvcs S «I? S sin (^) - V^t)V^) E S »^P sin ('^) | ( 28 ) 
^j2A^(e, = +l,P zz )sm(n^, (29) 

n=l 



A LL (e e =+l,P z 



1 



Ada Men = +l,P z 



da^+(P ZZl( f>)+da^+(P z 



da^+(P zz ,< 



da^+(P x 



A 



dcr TI g(e£ — +1, P z 



Pi (0)^2(0) 



X cf LP cos(t^) + A DVC s X c°I£ s cos(n0) 



n=0 



E^ 



cos(n4 
LL 



= +l,P zz )cos(n<) 



(30) 



(31) 
(32) 



The approximation in Eqs. 27, 29, and 32 is due to the truncation of terms in the 
Fourier series arising from the azimuthal dependences in the common denominator and 
the lepton propagators of Eqs. 26, 28, and 31. The Fourier coefficients of the expansion of 
the asymmetries are hereafter called asymmetry amplitudes. Although these asymmetry 
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amplitudes differ from the coefficients appearing in Eqs. 8-10 and Eqs. 26, 28, and 31, 
they may provide similar information in the comparison of model predictions with data. 



5.2. Single-beam-helicity formalism 

In order to extract more information on various combinations of Fourier coefficients in 
Eqs. 8-10, it is possible to use data collected with negative polarization of the e~ beam in 
conjunction with the subset of positron data with the same sign of the beam polarization. 
In this case, another set of Fourier coefficients of the single-beam-helicity asymmetries 
A c <=(Pe, Pzzi (f>); Ao (Pi,P XZ :<i>) an d Ac (Pe,P zz ,(f>) can be simultaneously extracted, 

C 

where the subscript <— => indicates the charge asymmetry for a lepton beam with negative 
polarization on a longitudinally polarized deuterium target with vanishing net vector 
polarization. The subscript <h-L indicates the asymmetry with respect to longitudinal 

vector target polarization for a charge-averaged lepton beam again with negative beam 

c 

polarization. Similarly, the subscript <— L indicates the double asymmetry with respect 
to lepton charge and longitudinal vector target polarization. 

The azimuthal distribution in the expectation value of the yield in this case can be 
written as 



d(A0(e/, P t , P z ,P Z zA) =C(et,Pe, P Z ,P ZZ ) v(<t>) da^(P e , P zz ,<f>) 

l + e e Ac^(Pe,P zz ,<f>) + PzAo T (Pz,Pzz,<t>) + e i P z Ac {P e , P zz ,<j>)} . (33) 



Here, the cross section da ^{Pi , P ZZ , 4>) for production of real photons by a charge- 
averaged polarized lepton beam on a tensor-polarized deuterium target with vanishing 
vector polarization is defined as 



daoAPe,P z 



1 



da^+ {Pi, P zz , <j>) + d<7^+ (Pt ,Pzz,< 



da^~(Pi, P zz ,0) + da^- {Pi, P zz , , 



K 



lBH 



2 ^ 2 

J2 C ™ C ° S (^) + P - 5>™ - Cn "LLP) C0S(7 



n=0 



K 



DVCS 



£ « P S cos(n0) + P e J2 -™ Mr 



n=0 
2 



\P ZZ ( E(Cn P S - cos(n0) + PiJ2 «™ S - S °X£p) M**) 



(34) 



(35) 



Then the single-beam-helicity asymmetries appearing in Eq. 33 are expressed as 
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x<^ da^+(P e ,P zz ,q>)+da^+(P e ,-P z ,P z 



da^~{P e , P zz ,^) + da^-(P e , P zz ,< 
K 



A, 



Vi(<i>)M<f>) 



da^(P e ,P zz ,(j)) 

3 2 

E c n,u„ P cos(n0) + P e 4,un P sin(n0) 



n=Q 



2 



9 P ^ ( H( C «,unP ~ C «,LLP) COs(n0) + P t ^( 

s n,unp — s n,LLp) sill(n</>) 



£ A C C °: ( "«(P^) cos(n0) + P, £ Afi" 0) (P zz ) sin(n0) 



n=0 



•4o,(P/,P* 



4dao $ (P^,P zz ,0) 

{ [da S +(P,,P ZZj( /») +da^-(P £ ,P z 

da^ + (P £ ,P 2Z ,0) +da^(P,,P Z2 , 
A" 



A 



BH 



da^{P e ,P zz ,4>) 
1 



n=Q 



A DVCS [P £ <I? S C ° S (-^) + E ^ Si *(^) 



n=0 



P t V A C OS( ^(P ZZ ) cos(n0) + V AT n *\P zz ) sin(n^) 



n=0 
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Ac (Pe,P zz ,< 



4dffog(P*,P„,0) 

| [da S +(P,,P zz ,</>) +da^-(P e ,P z 

dc^+(P £ ,P z ^) +da S ~(P,,P Z2 , 
K 



(42) 



LP cos(n0) 

2 



3 

S ",LP 

n=l 

3 



Pl(^)P 2 (0) 

2 

< — . — /__ i ■ 

(Pzz) cos(n 

< — L 

All the asymmetries defined in this paper are summarized in Table 3 



sin(n0) | 
^Af W) (P 2Z )sin(^) 



(43) 
(44) 



Table 3 

Extracted beam-helicity, beam-charge and target-spin asymmetries on a polarized deuterium target. 
The symbol ■ marks which data taken under certain experimental conditions (beam polarization, beam 
charge and target polarization state) are available for the construction of the respective asymmetry. The 
— or + indicates the sign with which the corresponding yield enters the numerator of the asymmetry. For 
the case that the target is populated with deuterons in the state A = ±1, the ideal target polarizations 
are P z = ±1 and P zz = 1, while for the case A = 0, P z = and P zz = —2. The sensitivity of coherent 
scattering to the corresponding Compton form factors or BH amplitude is indicated. 







Lepton charge 


Target population (deuterons) 


Beam helicity 










A = +1 A = -1 A = 


A = +1 A = — 1 


Coherent 






+1 -1 


=> -4= 


— > <- 


sensitivity 


CD 
aC 
H 


a^ 




■ + ■ 


■ ■ 


$Sm(Hi,H B ) 


S ingle- cha 


Ajl 






■ + ■ 


Qm(Hi) 


All 






■ - ■ 


(BH) 






■ + ■ ■ 


■ - ■ 




>. 
"0 




■ - ■ 


■ + ■ 




Sm/SRe(«i,H 5 ) 


"o 


A 


■ + ■ 






(BH) 


3ingL 


A c 

< — L 


■ - ■ 






9ta/SRe(«i) 



6. Background corrections and systematic uncertainties 

The asymmetry amplitudes are corrected for background contributions, mainly decays 
to two photons of semi-inclusive neutral mesons, using the method described in detail in 
Ref. [42]. The average contribution from semi-inclusive background is 4.6%. The contri- 
bution of exclusive pions is neglected, as it is found to be less than 0.7% in each kinematic 
bin, supported by studies of HERMES data [43] . After applying this correction, the re- 
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suiting asymmetry amplitudes are expected to originate from coherent and incoherent 
photon production, the latter possibly including nucleon excitation. 

The dominant contributions to the total systematic uncertainty are the effects of the 
limited spectrometer acceptance and from the finite bin widths used for the final presen- 
tation of the results. The latter originates from the difference of the amplitudes integrated 
over one bin in all kinematic variables, compared to the asymmetry amplitudes calcu- 
lated at the average values of the kinematic variables. The combined contribution to the 
systematic uncertainty from limited spectrometer acceptance, finite bin width, and the 
alignment of the spectrometer elements with respect to the beam is determined from a 
Monte Carlo simulation using a convenient parameterization [44] of the VGG model [45] 
(sec details in Ref. [27]). Five GPD model variants are considered, including only incoher- 
ent processes on the proton and neutron. In each kinematic bin, the resulting systematic 
uncertainty is defined as the root-mean-square average of the five differences between the 
asymmetry amplitude extracted from the Monte Carlo data and the corresponding model 
predictions calculated analytically at the mean kinematic values of that bin. In the case 
of the single-charge beam-helicity asymmetry, all five models overpredict the magnitudes 
of the sin <fi harmonics by about a factor of two, leading to a probable overestimate of this 
contribution to the uncertainties. The other source of uncertainty is associated with the 
background correction. For asymmetries involving target vector polarization, no system- 
atic uncertainty due to luminosity is assigned. This is legitimate because the luminosity 
does not depend on the target polarization, the target polarization flips rapidly compared 
to changes in luminosity, and beam polarization dependent weights are assigned to each 
event in the extraction. There is an additional overall scale uncertainty arising from the 
uncertainty in the measurement of the beam and/or target polarizations. Not included is 
any contribution due to additional QED vertices, as for the case of polarized target and 
polarized beam the most significant of these has been estimated to be negligible [46] . The 
total systematic uncertainty in a kinematic bin is determined by adding quadratically all 
contributions to the systematic uncertainty for that bin. 

7. Results 

7.1. Single- and double-spin asymmetries 

The results for the Fourier amplitudes of the single-charge asymmetries A^(eg = 
+l,P zz ,4>), Auhiee = +l,P zz ,4>) and A^(ei> = +1,P ZZ ,4>) are presented in Figs. 1- 
3 as a function of — t, xn, or Q 2 and are also given in Table 5. While the variable 
Xu would be the appropriate choice to present experimental results for pure coherent 
scattering, the nucleonic Bjorken variable a; at is the practical choice in this case where 
incoherent scattering dominates over most of the kinematic range. The 'overall' results 
in the left columns correspond to the entire HERMES kinematic acceptance. Figure 1 
shows the amplitudes A™^ n ^ related to beam helicity only, while Figs. 2 and 3 show the 
amplitudes A*^ n ^ , which relate to target vector polarization only, and the amplitudes 
^■ll ' wn ich relate to the product of beam helicity and target vector polarization. 
Table 4 and Fig. 4 show in each kinematic bin the estimated fractional contributions 
to the yield from the coherent process and from processes leading to baryonic resonant 
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Fig. 1. Results from the present work (red filled squares) representing single-charge beam-helicity asym- 
metry amplitudes X ai ^ ' describing the dependence of the sum of squared DVCS and interference terms 
on the beam helicity, for a tensor polarization of P zz = 0.827 (indicated by the symbol <->). The black 
open squares represent charge-difference amplitudes A^y"^ from only the interference term, extracted 
from unpolarized deuterium data [27]. The error bars represent the statistical uncertainties, while the 
coarsely hatched (open) bands represent the systematic uncertainties of the filled (open) squares. There 
is an additional overall 1.9% (2.4%) scale uncertainty arising from the uncertainty in the measurement of 
the beam polarization in the case of polarized (unpolarized) deuterium data. The points for unpolarized 
deuterium data are slightly shifted to the left for better visibility. The finely hatched band shows the 
results of theoretical calculations for the combination of incoherent scattering on proton and neutron, 
using variants of the VGG double-distribution model [45,48] with a Regge ansatz for modeling the t 
dependence of GPDs [49] . 



final states. They are obtained from a Monte Carlo simulation using an exclusive-photon 
generator described in Ref. [27]. 

The values for the sin <p amplitude of the asymmetry A^ in Fig. 1 are found to be 
significantly negative, while the sin(20) amplitude is found to be consistent with zero. 
Figure 1 also presents for comparison the amplitudes of the charge-difference asymmetry 
A\^j extracted from a previous measurement on unpolarized deuterons [27]. Under the 
same approximations as those leading to Eq. 13, A^% is expected to differ from A\-^ 
(only if P zz 0) due only to a term involving the CFF H^. Figure 1 shows that these 
two asymmetries are found to be consistent in most kinematic regions, except possibly 
for the last —t or xn bin in the case of sin(2</>). (The overall results differ by only 1.7 
standard deviations in the total experimental uncertainties. ) The consistency in the first 



Here and hereafter we neglect any possible correlations arising from common treatments of different 
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-0.4 ■ — ; ■ ■ ^^^^^^^^ = 

10" 1 10" 1 1 10 

overall -t [GeV 2 ] x N Q 2 [GeV 2 ] 

Fig. 2. Single-charge target-spin asymmetry amplitudes describing the dependence of the sum of squared 
DVCS and interference terms on the target vector polarization, for a tensor polarization of P zz = 0.827. 
The squares represent the results from the present work. The triangles denote the corresponding am- 
plitudes extracted from longitudinally polarized hydrogen data [24], The error bars (bands) represent 
the statistical (systematic) uncertainties. The finely hatched bands have the same meaning as in Fig. 1. 
There is an additional overall 4.0% (4.2%) scale uncertainty arising from the uncertainty in the mea- 
surement of the target polarization in the case of deuterium (hydrogen). The points for hydrogen are 
slightly shifted to the left for better visibility. 

—t bin, where the contribution from coherent scattering is significant, suggests that there 
is no distinctive contribution from %s , as was observed in the case of the corresponding 
forward limit [20,21]. 

In the first —t bin, the asymmetry amplitude ^""^oh ^ or P ure coherent scattering on 
a polarized deuterium target was estimated from the measured asymmetry by correcting 
for the incoherent contributions of the proton and neutron and their resonances (sec 
Ref. [27]). This correction is based on the assumption that for the incoherent contribution 
of the proton, A S ^^(P ZZ = 0.827) « ^lu i wnere the latter was measured on a hydrogen 
target [26]. The fractional contributions and the asymmetry for incoherent scattering 
from the neutron was taken from the Monte Carlo calculation described in section 6, 
with uncertainties equal to their magnitude. The result for the asymmetry amplitude 
A tmtoh( P z* = °- 827 ) is estimated t0 be -°- 12 ± 0.17(stat.) ± 0.14(syst.) ± 0.02(model), 
where the systematic uncertainty is propagated from only the corresponding experimental 
uncertainties. Within the uncertainties there is no evidence of a difference between this 
value and the value for the asymmetry amplitude -AluI coh = — 0.29 ± 0.18(stat.) ± 
0.03(syst.) previously estimated for coherent scattering on an unpolarized deuterium 



data sets. 
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Fig. 3. Single-charge double-spin asymmetry amplitudes describing the dependence of the sum of 
Bcthe-Heitlcr, squared DVCS and interference terms on the product of the beam helicity and target 
vector polarization, for a tensor polarization of P zz = 0.827. The plotted symbols and bands have the 
same meaning as in Fig. 2. There is an additional overall 4.4% (5.3%) scale uncertainty arising from 
the uncertainties in the measurement of the beam and target polarizations in the case of deuterium 
(hydrogen) data. 
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Fig. 4. Simulated yield fractions of coherent and resonant production 

target, using a disjoint HERMES data set for an unpolarizcd deuterium target [27], but 
using the same data set for a hydrogen target. 

The extracted values for the sine/) and sin(2^) amplitudes of the single-charge asym- 
metry ^4ul measured on a longitudinally polarized deuterium target are shown in Fig. 2. 
The 'overall' values are slightly negative by less than 1.5 standard deviations of the total 
experimental uncertainty. For coherent scattering on the deuteron, the amplitude ^ul^ 
is sensitive to the imaginary part of a combination of deuteron CFFs Hi, Hi and 
weighted with the elastic form factors of the deuteron G\ and G2 (see Eq. 21). In particu- 
lar, for the first —t bin where (xd) = 0.04, G\ is about 30 times larger than ^rr-Gz- Thus 
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Table 4 

Simulated fractional contributions of coherent and resonant processes on a deuteron, in each kinematic 
bin. 
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the CFF T-Li may influence the resulting A s ^ amplitude in the first —t bin where the 
coherent process contributes approximately 40%. For comparison, the same amplitudes 
measured on a longitudinally polarized hydrogen target [24] are also shown in Fig. 2. The 
sin <fi amplitude shows consistency between deuterium and hydrogen data both for the 
'overall' result and the kinematic projections on — t, xn, and Q 2 . In this comparison, no 
account was taken of the 7.5% depolarization of nucleons in the deuteron due to the 5% 
admixture of the D-state [47]. The 'overall' results on the sin(20) amplitude differ be- 
tween the two targets by 1.5 standard deviations of the total experimental uncertainties, 
mainly due to the region of large —t, but in only one xn bin. The 'overall' result on the 
asymmetry amplitude A^j 3 ^ is slightly negative by less than 1.7 standard deviations 
of the total experimental uncertainty. The sin(3</)) amplitude shows consistency between 
deuterium and hydrogen data, accounting for the total experimental uncertainties of the 
corresponding measurements, except possibly for the highest xn bin. 

The ^4^( n< ^ amplitudes of the single-charge double-spin asymmetry measured using 
longitudinally polarized deuteron data and presented in Fig. 3 are found to be compatible 
with zero, although the amplitude is positive by 1.6 standard deviations of the total 

experimental uncertainty. Within the uncertainties, these asymmetry amplitudes do not 
show significant differences from those measured on a longitudinally polarized hydrogen 
target [24] , except possibly for the overall result for the amplitude A™^ ^ ; where there is 
observed a discrepancy of 1.9 standard deviations in the total experimental uncertainties. 

The finely hatched bands in Figs. 1-3 represent results of theoretical calculations based 
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on the GPD model described in Ref. [45] , using the VGG computer program of Rcf. [48] . 
The Regge ansatz for modeling the t dependence of GPDs [49] is used in these calcu- 
lations. The model [45] is an implementation of the double-distribution concept [1,2] 
where the kernel of the double distribution contains a profile function that determines 
the dependence on £, controlled by a parameter b [50] for each quark flavor. The cross 
sections are calculated as the sum of the incoherent processes on the proton and neutron 
in each kinematic bin. (No computer program is available simulating coherent scattering 
on the deuteron.) The width of the theoretical bands in Figs. 1-3 corresponds to the 
range of values of the asymmetry amplitudes obtained by varying the profile parameters 
& V ai and & sea between unity and infinity. In the comparison of these predictions with ex- 
perimental results, it should be noted that the effect of the D-state of the deuteron on 
the polarization of the nuclcons inside the deuteron was not taken into account. 

The model calculations predict a magnitude of the sin <p harmonic of the single-charge 
beam-helicity asymmetry that exceeds that of the data by about a factor of two, a 
situation similar to that found in the case of a hydrogen target [26] . On the other hand the 
predictions arc in good agreement with data for single-charge target-spin asymmetries. A 
large difference appears between the predictions for the sin <f> harmonic of this asymmetry 
on the deuteron and proton targets, arising entirely from the contributions of the neutron. 
The data are consistent with this difference, but lack the precision to confirm the large 
positive prediction of the neutron asymmetry by this model. The predictions are in 
good agreement with the single-charge double-spin asymmetry amplitudes, aside from 
the cos(0</>) harmonic. Here the theoretical predictions for both the deuteron and proton, 
which are dominated by the BH contribution, are significantly positive, in agreement 
with the proton data, while the more precise deuteron data are consistent with zero. The 
small contribution of coherent scattering to the overall result, with a predicted negative 
asymmetry [18], is expected to slightly reduce this asymmetry amplitude for the deuteron. 

7.2. The beam-charge, charge-averaged, and beam- charge® target- spin asymmetries 

The results for the Fourier amplitudes of the single-bcam-hclicity asymmetries are 
presented in Figs. 5-7. More specifically, Figs. 5, 6, and 7 show the cos(n</>) and sin(n0) 
harmonics of the asymmetry Ac<=(Pe, Pzz,4>)> Ao (Pe,P zz ,4>) and Ac (Pe, P zz , 4>), re- 
spectively (see also Tables 6-8)3^ P e = -0.530 ± L 0.012 and P zz = 0.827 ± 0.027. 

The only overall results for the asymmetry Ac^ in Fig. 5 that are found to be sig- 
nificantly non-zero are the cos cj) and sin <f> amplitudes. The theoretical calculations for 
incoherent scattering predict that the results for the amplitudes A c ° s J; n ^ should strongly 

resemble those for the amplitudes A c £ s<Jl ^ measured with an unpolarized beam on an 
unpolarized deuterium target [27]. The data confirm this resemblance, even in the first 
—t bin where coherent scattering contributes about 40% of the yield. This is another 
indication that the CFF [18], in this case its real part, makes no distinctive contri- 
bution to coherent scattering off deuterons, similar to the case of , as was noted in 
the discussion in Section 7.1 about the dependence of of Eq. 13 on the imaginary 

part of this CFF. 

The numerators of the A*™^ 1 ^ amplitudes shown in Fig. 5 differ from those of the 
sin(n<f>) amplitudes of the A^ asymmetry shown in Fig. 1 only by squared DVCS terms. 
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Fig. 5. Results from the present work (red filled squares) representing single-beam-helicity charge asym- 
metry amplitudes J 4**j( n 'M ant j for P e = —0.530 and a tensor polarization of P zz = 0.827 

(indicated by the symbol The black open squares are J 4^ os (™0) amplitudes extracted from data 
recorded with an unpolarized beam and unpolarized deuterium target [27]. The error bars and bands 
and finely hatched bands have the same meaning as in Fig. 1 . The points for unpolarized deuterium data 
are slightly shifted to the left for better visibility. There is an additional overall 2.2% scale uncertainty 
for the A 8 ^? amplitudes arising from the uncertainty in the measurement of the beam polarization. 
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Furthermore, the cross sections da <= and dcr^ in the denominators of these two asym- 
metries should be similar because they are dominated by Bethc-Hcitlcr contributions. 
Hence, these asymmetry amplitudes are expected to be similar, and within the statisti- 
cal accuracy this is indeed found to be the case. 
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Fig. 6. Kinematic dependence of the charge-averaged single-beam-hclicity target-spin asymmetry am- 
plitudes A c o os(n</>) and A 8 ^"^, for P e = -0.530 and a tensor polarization of P xx = 0.827. The plotted 

symbols and bands have the same meaning as in Fig. 5. There is an additional overall 5.3% (5.7%) 

( A c ° s ^ n ^j amplitudes arising from the uncertainties in the 



scale uncertainty for the extracted A a ^ n 



measurement of the target (beam and target) polarizations. 



The cos(n</>) amplitudes of the asymmetry Ao in Fig. 6 contain a sum of BH and 
squared DVCS even harmonics, and relate to the longitudinal vector polarization of 
the target. However, even where the BH contribution dominates the numerator of the 
asymmetry amplitude A c ° s< " ^ for incoherent scattering at not small —t, the data are 
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found to be consistent with zero, and differing by 1.7 standard deviations in the total 
experimental uncertainty from the positive prediction for the overall result. The sin(n0) 
amplitudes of the asymmetry A o in Fig. 6 receive contributions from the pure squared 
DVCS harmonics only, and arc found to be consistent with zero. 

Of particular interest are the A c ° s ( n ^ a nd ^ n ( ni W amplitudes shown in Fig. 7, which 

represent respectively the even and odd vector-polarization related harmonics of the 
interference term only, receiving no contribution from pure BH and DVCS terms. The- 
theoretical predictions for the cos(n^>) harmonics are negligibly small, while the data differ 
from zero by about two standard deviations for the first two harmonics. As expected and 
observed in the case of unpolarized hydrogen and deuterium targets, the cos(0(f>) and 
cos 4> harmonics are found to have opposite signs. 

Like the asymmetry amplitude , m the first —t bin the asymmetry amplitude 

is sensitive to the imaginary part of the deuteron CFF Hi- Within their statistical 

accuracies, they are found to be consistent, although receives also a contribution 

from the squared DVCS term (see Eq. 29). The asymmetry amplitude ^4™ b * is sensitive 

to the real part of the deuteron CFF Hi- Unlike the corresponding harmonic ^li!^' ^ 
does not receive a contribution from the Bethe-Heitler term. The sin(n0) harmonics are 
found to be consistent with zero and also with the small negative prediction in the case 
of the sin</> harmonic. 

From the definitions of the asymmetries -4ul, All, ■Ao and Ac in Eqs. 19, 30, 39, 
and 42, and also from examination of Table 2, it can be seen that they are related. In the 
case of approximate equality of dtr ^ and da v ^, the following relations hold between 
the asymmetry amplitudes: 

1,2, (45) 
0, 1 . (46) 

For most of the kinematic points, the differences between left and right hand sides of Eqs. 
45 and 46 are found below 1.2 standard deviations of the total experimental uncertainties, 
while for the remaining six points they are between 1.5 and 2.0. Note that here the 
correlations between two asymmetries from the right hand sides are taken into account. 

7.3. The beam-helicity <g> tensor asymmetry Alzz 

The definition of the asymmetry Alzz is given in Eq. 16. As mentioned in Section 3, 
for the extraction of this asymmetry, the data taken with a positron beam and with 
the average target tensor polarization P zz = —1.656 are used in combination with the 
positron data collected on a longitudinally polarized deuterium target with P zz = 0.827. 
The same maximum likelihood technique [42] unbinned in azimuthal angle c/> was used 
to extract the .A^i Fourier amplitudes. The A^f amplitude is found to have the 
value .A L ™ = — 0.130 ± 0.121(siai.) ± 0.051(sj/s<.) when extracted in the entire kine- 
matic range of the data set, while in the region —t < 0.06 GeV 2 , where the contribution 
from coherent scattering on a longitudinally polarized deuteron is approximately 40%, 
this value is found to be 0.074 ± 0.196(stat.) ± 0.022(syst.). These results are subject 
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^UL — ^0 t A C l n — 
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Fig. 7. Kinematic dependence of the single-beam-helicity beam-chargeCS>target-spin asymmetry ampli- 



tudes A°„ o8(n0) and AT (r ' 



for 



-0.530 and a tensor polarization of P zz = 0.827. The plotted 



symbols and bands have the same meaning as in Fig. 5. There is an additional overall 5.3% (5.7%) 
scale uncertainty for the extracted j{ s ^ n ^ ^cos(n^)j amplitudes arising from the uncertainties in the 
measurement of the target (beam and target) polarizations. 

to an additional scale uncertainty of 3.8% arising from beam and target-tensor polariza- 
tions. The Fourier amplitudes related to higher twist are found to be compatible with 
zero within the statistical uncertainties. This 'zero' result for the beam-helicity<g>tensor 
asymmetry Alzz extracted independently of the results for A^u and *4 L g confirms that 
there is no distinctive contribution from the dcuteron CFF H 5 for coherent scattering. 



8. Summary 

Azimuthal asymmetries with respect to target polarization alone and also combined 
with beam hclicity and/or beam charge for hard exclusive electroproduction of real pho- 
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tons in decp-inclastic scattering from a longitudinally polarized deuterium target are 
measured for the first time. The asymmetries are attributed to the interference between 
the deeply virtual Compton scattering and Bcthe Hcitlcr processes. The asymmetries arc 
observed in the exclusive region -(1.5) 2 GcV 2 < Mf < (1.7) 2 GeV 2 of the squared miss- 
ing mass. The dependences of these asymmetries on —t, xjv, or Q 2 are investigated. The 
results include the coherent process ed ~ > ec?7 and the incoherent process ed-> epnj 
where in addition a nucleon may be excited to a resonance. Within the total experimental 
uncertainties, the results of the sinusoidal (cosinusoidal) amplitudes of the asymmetry 
v4 L ^ (A c ^) extracted from a data set with P zz = 0.827 (corresponding to a small popu- 
lation for the A = state) resemble those for the amplitudes extracted from unpolarizcd 
deuterium data at HERMES. Therefore, no indication of effects of tensor polarization 
was found at small values of —t, in particular in the region — t < 0.06 GeV 2 where the 
coherent process contributes up to 40%. Neither the A^j nor amplitudes 
measured on longitudinally polarized deuterons show significant differences compared 
with those extracted from longitudinally polarized protons, considering the total exper- 
imental uncertainties. (Statistically marginal differences are observed for A^} 2 ^ and 

^cos(O0)-j 

The sinusoidal amplitudes of the tensor asymmetry Al Z z are compatible with zero for 
the whole kinematic range as well as for the region —t < 0.06 GeV 2 within the accuracy 
of the measurement. This suggests that differences between the leading amplitudes of 
the asymmetries -4ltj anc ^ A^ for coherent scattering from unpolarizcd and longitudi- 
nally polarized deuterium targets, respectively, should be small. Indeed, within the total 
experimental uncertainties, no difference is seen between the reconstructed values of the 
asymmetry amplitudes ^"toh anc ^ ^Lut coh- 

In conclusion, even in the region —t < 0.06 GeV 2 where the coherent process con- 
tributes about 40%, all asymmetries on deuterium that have (approximate) counterparts 
for hydrogen are found to be compatible with them. The data are unable to reveal any 
evidence of the influence of the Compton form factor H5 or features of the deuteron 
Compton form factors Hi and Hi that distinguish them from the counterparts for the 
proton. Hence, coherent scattering presents no obvious signature in these data. The 
deuteron Compton form factor Hi appears to have a similar behavior as H of the pro- 
ton. The data were compared with theoretical calculations for only incoherent scattering, 
based on a well-known GPD model. Those asymmetries that are expected to resemble 
counterparts for a hydrogen target reveal the same shortcomings of the model calculations 
that appeared in comparisons with the hydrogen data. 
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Table 5 

Results for azimuthal Fourier amplitudes of the single-charge asymmetries -4 L ^, ytijL and All, extracted 
from longitudinally polarized deuteron data, for a tensor polarization of P zz = 0.827. Not included are 
the 1.9%, 4.0% and 4.4% scale uncertainties for corresponding asymmetry amplitudes arising from the 
uncertainties in the measurement of the beam, target, beam and target polarizations, respectively. 
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Table 6 

Results for azimuthal Fourier amplitudes of the single-beam- helicity charge asymmety >t c<= , extracted 
from longitudinally polarized deuteron data, for P( = —0.530 and a tensor polarization of P zz = 0.827. 
Not included is the 2.2% scale uncertainty for sin(n0) asymmetry amplitudes arising from the uncertainty 
in the measurement of the beam polarization. 



Kinematic bin <-t> (x N ) (Q 2 ) A?'"*' A"°"J A ^(^) 

[GcV 2 ] [GcV 2 ] ±5 sta t ± <5 5yBt ±«rtat ± S eyat ±«st»t ± <W=t 
Overall 0.13 0.10 2.5 -0.012 ± 0.018 ± 0.034 0.065 ± 0.026 ± 0.009 0.017 ± 0.026 ± 0.003 
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0.00-0.06 0.03 0.08 1.9 0.006 ± 0.030 ± 0.031 0.001 ± 0.041 ± 0.012 0.005 ± 0.042 ± 0.009 
0.06-0.14 0.10 0.10 2.5 0.074 ± 0.035 ± 0.034 0.037 ± 0.049 ± 0.008 0.046 ± 0.049 ± 0.006 
0.14-0.30 0.20 0.11 2.9 -0.098 ± 0.036 ± 0.031 0.139 ± 0.052 ± 0.008 -0.007 ± 0.051 ± 0.005 
0.30-0.70 0.42 0.12 3.5 -0.086 ± 0.058 ± 0.028 0.159 ± 0.088 ± 0.007 0.056 ± 0.079 ± 0.008 
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0.03-0.07 0.11 0.05 1.4 -0.046 ± 0.031 ± 0.035 0.042 ± 0.044 ± 0.009 -0.078 ± 0.043 ± 0.005 
0.07-0.10 0.11 0.08 2.1 0.025 ± 0.035 ± 0.030 0.005 ± 0.050 ± 0.005 0.071 ± 0.048 ± 0.003 
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> 



T 


0.00-0.06 0.03 0.08 1.9 -0.018 ± 0.042 ± 0.004 -0.158 ± 0.081 ± 0.050 0.109 ± 0.080 ± 0.005 
0.06-0.14 0.10 0.10 2.5 0.075 ± 0.049 ± 0.004 -0.156 ± 0.095 ± 0.036 -0.021 ± 0.094 ± 0.015 
0.14-0.30 0.20 0.11 2.9 0.053 ± 0.052 ± 0.005 -0.126 ± 0.098 ± 0.021 0.018 ± 0.099 ± 0.011 
0.30-0.70 0.42 0.12 3.5 0.098 ± 0.077 ± 0.007 0.141 ± 0.142 ± 0.029 -0.015 ± 0.153 ± 0.015 




0.03-0.07 0.11 0.05 1.4 -0.002 ± 0.041 ± 0.002 -0.109 ± 0.080 ± 0.060 -0.028 ± 0.081 ± 0.039 
0.07-0.10 0.11 0.08 2.1 0.028 ± 0.049 ± 0.003 -0.069 ± 0.092 ± 0.059 0.095 ± 0.094 ± 0.018 
0.10-0.15 0.14 0.12 3.1 0.091 ± 0.056 ± 0.006 -0.288 ± 0.107 ± 0.058 0. 119 ± 0. 107 ± 0.018 
0.15-0.35 0.20 0.20 5.0 0.089 ± 0.075 ± 0.005 -0.032 ± 0.147 ± 0.054 -0.020 ± 0.141 ± 0.011 
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1.0-1.5 0.09 0.06 1.2 -0.020 ± 0.048 ± 0.002 0.020 ± 0.095 ± 0.041 0.011 ± 0.092 ± 0.056 
1.5-2.3 0.11 0.08 1.9 0.074 ± 0.047 ± 0.003 -0.266 ± 0.087 ± 0.043 0.095 ± 0.092 ± 0.038 
2.3-3.5 0.14 0.11 2.8 0.076 ± 0.053 ± 0.006 -0.076 ± 0.101 ± 0.053 -0.023 ± 0.099 ± 0.016 
3.5-10.0 0.20 0.17 4.9 0.059 ± 0.058 ± 0.003 -0.145 ± 0.112 ± 0.044 0.042 ± 0. 112 ± 0.011 
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Table 7 

Results for azimuthal Fourier amplitudes of the singlc-bcam-hclicity charge-averaged asymmetry A o , 
extracted from longitudinally polarized deuteron data, for P( = —0.530 and a tensor polarization of P zz = 
0.827. Not included is the 5.3% (5.7%) scale uncertainty for the sva(n<j>) (cos(n<j>)) asymmetry amplitudes 
arising from the uncertainties in the measurement of the target (beam and target) polarizations. 



Kinematic bin 


(-*) 
[GoV 2 ] 




(Q 2 > 

[GoV 2 ] 


COB (0*) 


±<5 s tat i <5syst 


«— L 

=t"5 B tat ± <5 ays t 


.sin (2*) 




Overall 


0.13 


0.10 


2.5 


0.021 ± 0.044 ± 0.009 


-0.041 ± 0.062 ± 0.010 


0.005 ± 0.033 ± 0.003 


-0.036 ± 0.033 ± 0.003 




0.00-0.06 


0.03 


0.08 


1.!) 


-0.009 ± 0.072 ± 0.008 


0.087 ± 0.101 ± 0.011 


0.011 ± 0.054 ± 0.004 


0.033 ± 0.054 ± 0.005 


> 


T 


0.06-0.14 


0.10 


0.10 


2.5 


0.039 ± 0.083 ± 0.012 


-0.005 ± 0.115 ± 0.011 


0.003 ± 0.063 ± 0.006 


-0.108 ± 0.062 ± 0.006 


0.14-0.30 


0.20 


0.11 


2.!) 


0.030 ± 0.091 ± 0.008 


-0.282 ± 0.128 ± 0.024 


0.058 ± 0.068 ± 0.003 


-0.034 ± 0.068 ± 0.005 




0.30-0.70 


0.42 


0.12 


3.5 


0.024 ± 0.142 ± 0.014 


-0.056 ± 0.206 ± 0.059 


-0.137 ± 0.100 ± 0.008 


-0.093 ± 0.111 ± 0.011 




0.03-0.07 


0.11 


0.05 


1.4 


0.051 ± 0.075 ± 0.004 


-0.121 ± 0.106 ± 0.005 


-0.001 ± 0.054 ± 0.002 


-0.008 ± 0.056 ± 0.002 


z 


0.07-0.10 


0.11 


0.08 


2.1 


-0.014 ± 0.083 ± 0.012 


-0.049 ± 0.117 ± 0.021 


0.020 ± 0.062 ± 0.004 


-0.043 ± 0.062 ± 0.006 




0.10-0.15 


0.14 


0.12 


3.1 


-0.158 ± 0.095 ± 0.019 


0.033 ± 0.133 ± 0.030 


0.093 ± 0.072 ± 0.006 


-0.029 ± 0.071 ± 0.006 




0.15-0.35 


0.20 


0.20 


5,0 


0.228 ± 0.135 ± 0.027 


0.093 ± 0.187 ± 0.044 


-0.108 ± 0.100 ± 0.011 


-0.069 ± 0.095 ± 0.011 




1.0-1.5 


0.09 


0.06 


1.2 


0.028 ± 0.085 ± 0.009 


-0.108 ± 0.116 ± 0.013 


-0.011 ± 0.065 ± 0.008 


0.019 ± 0.065 ± 0.005 


> 



a 


1.5-2.3 


0.11 


0.08 


1.0 


-0.029 ± 0.082 ± 0.006 


-0.077 ± 0.119 ± 0.009 


-0.025 ± 0.059 ± 0.005 


-0.140 ± 0.061 ± 0.005 


2.3-3.5 


0.14 


0.11 


2.8 


0.018 ± 0.091 ± 0.012 


-0.049 ± 0.125 ± 0.013 


0.047 ± 0.069 ± 0.004 


-0.011 ± 0.066 ± 0.005 


Qf 


3.5-10.0 


0.20 


0.17 


4.9 


0.078 ± 0.100 ± 0.022 


0.127 ± 0.142 ± 0.020 


0.024 ± 0.075 ± 0.005 


-0.013 ± 0.075 ± 0.007 
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Table 8 

Results for azimuthal Fourier amplitudes of the single-beam-hclicity beam-chargc<X>targct-spin asym- 
metry A c , extracted from longitudinally polarized deuteron data, for Pi = —0.530 and a tensor po- 

4— L 

larization of P zz = 0.827. Not included is the 5.3% (5.7%) scale uncertainty for the sin(n</)) (cos(n<f>)) 
asymmetry amplitudes arising from the uncertainties in the measurement of the target (beam and target) 
polarizations. 



Kinematic bin 


<-*> 
[GeV 2 ] 


(*n) 


(Q 2 ) 

[GcV 2 ] 


.cos (04.) 

±<5 s tat i <5 s yst 


. cos cj> 

±5 a tat i ^syst 


.cos (20) 

±<5stat ± <5 s yst 




Overall 


0.13 


0.10 


2.5 


-0.082 ± 0.044 ± 0.002 


0.148 ± 0.062 ± 0.007 


-0.044 ± 0.057 ± 0.002 




0.00-0.06 


0.03 


0.08 


1.9 


-0.089 ± 0.072 ± 0.002 


0.124 ± 0.100 ± 0.006 


-0.161 ± 0.095 ± 0.003 


> 



T 


0.06-0.14 


0.10 


0.10 


2.5 


-0.071 ± 0.082 ± 0.002 


0.147 ± 0.114 ± 0.007 


-0.225 ± 0.110 ± 0.006 


0.14-0.30 


0.20 


0.11 


2.9 


-0.152 ± 0.091 ± 0.003 


0.262 ± 0.127 ± 0.008 


0.092 ± 0.114 ± 0.005 




0.30-0.70 


0.42 


0.12 


3.5 


0.190 ± 0.142 ± 0.014 


0.244 ± 0.205 ± 0.006 


0.547 ± 0.172 ± 0.014 




0.03-0.07 


0.11 


0.05 


1.4 


-0.044 ± 0.075 ± 0.003 


0.101 ± 0.109 ± 0.001 


0.019 ± 0.093 ± 0.003 




0.07-0.10 


0.11 


0.08 


2.1 


-0.170 ± 0.082 ± 0.004 


0.210 ± 0.116 ± 0.006 


-0.132 ± 0.108 ± 0.003 




0.10-0.15 


0.14 


0.12 


3.1 


0.058 ± 0.095 ± 0.002 


0.323 ± 0.132 ± 0.010 


-0.239 ± 0.123 ± 0.007 




0.15-0.35 


0.20 


0.20 


5.0 


-0.281 ± 0.135 ± 0.012 


-0.045 ± 0.188 ± 0.009 


0.269 ± 0.173 ± 0.012 




1.0-1.5 


0.09 


0.06 


1.2 


-0.100 ± 0.085 ± 0.002 


-0.016 ± 0.116 ± 0.001 


-0.032 ± 0.108 ± 0.005 


> 

o 

o 


1.5-2.3 


0.11 


0.08 


1.9 


-0.097 ± 0.082 ± 0.002 


0.140 ± 0.119 ± 0.005 


-0.069 ± 0.104 ± 0.003 


2.3-3.5 


0.14 


0.11 


2.8 


-0.002 ± 0.089 ± 0.001 


0.483 ± 0.123 ± 0.010 


-0.068 ± 0.120 ± 0.002 


0- 


3.5-10.0 


0.20 


0.17 


4.9 


-0.145 ± 0.100 ± 0.003 


-0.003 ± 0.141 ± 0.008 


0.005 ± 0.131 ± 0.009 




Kinematic bin 


(-*> 
[GcV 2 ] 




(Q 2 > 

[GeV 2 ] 


=t<5 a tat i <5 sys t 


.sin (20) 

±<5stat i <5 a yst 


.sin (30) 

±<5 s tat i ^syst 




Overall 


0.13 


0.10 


2.5 


-0.023 ± 0.033 ± 0.028 


-0.035 ± 0.033 ± 0.008 


-0.009 ± 0.030 ± 0.003 




0.00-0.06 


0.03 


0.08 


1.9 


-0.032 ± 0.054 ± 0.033 


-0.116 ± 0.053 ± 0.007 


-0.064 ± 0.050 ± 0.004 


> 

a 


0.06-0.14 


0.10 


0.10 


2.5 


0.016 ± 0.062 ± 0.016 


-0.016 ± 0.062 ± 0.009 


-0.033 ± 0.059 ± 0.005 


0.14-0.30 


0.20 


0.11 


2.9 


-0.045 ± 0.068 ± 0.010 


0.024 ± 0.067 ± 0.015 


0.025 ± 0.060 ± 0.008 


i 


0.30-0.70 


0.42 


0.12 


3.5 


-0.001 ± 0.102 ± 0.005 


0.212 ± 0.115 ± 0.020 


0.201 ± 0.099 ± 0.008 




0.03-0.07 


0.11 


0.05 


1.4 


-0.027 ± 0.054 ± 0.015 


0.006 ± 0.058 ± 0.005 


-0.007 ± 0.051 ± 0.001 




0.07-0.10 


0.11 


0.08 


2.1 


-0.073 ± 0.061 ± 0.023 


-0.098 ± 0.062 ± 0.004 


0.001 ± 0.057 ± 0.003 




0.10-0.15 


0.14 


0.12 


3.1 


0.031 ± 0.072 ± 0.026 


-0.087 ± 0.072 ± 0.003 


0.025 ± 0.065 ± 0.003 




0.15-0.35 


0.20 


0.20 


5.0 


0.009 ± 0.100 ± 0.031 


-0.001 ± 0.097 ± 0.002 


-0.036 ± 0.090 ± 0.006 



1.0 


1.5 


0.09 


0.06 


1.2 


-0.003 ± 0.065 ± 0.022 


0.020 ± 0.065 ± 0.003 


-0.003 ± 0.058 ± 0.001 


1.5 


2.3 


0.11 


0.08 


1.9 


-0.064 ± 0.059 ± 0.026 


-0.041 ± 0.061 ± 0.004 


-0.017 ± 0.057 ± 0.005 


2.3 


3.5 


0.14 


0.11 


2.8 


-0.087 ± 0.068 ± 0.021 


-0.090 ± 0.065 ± 0.008 


-0.014 ± 0.062 ± 0.001 


3.5- 


10.0 


0.20 


0.17 


4.9 


0.102 ± 0.075 ± 0.027 


-0.040 ± 0.075 ± 0.005 


-0.004 ± 0.070 ± 0.004 
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